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The main hurdle in the realization of dynamic spectrum access (DSA) systems from physical 
layer perspective is the reliable sensing of low power licensed users. One such scenario shows up 
in the unlicensed use of TV bands where the TV Band Devices (TVBDs) are required to sense 
extremely low power wireless microphones (WMs). The lack of technical standard among various 
wireless manufacturers and the resemblance of certain WM signals to narrow-band interference signals, 
such as spurious emissions, further aggravate the problem. Due to these uncertainties, it is extremely 
difficult to abstract the features of WM signals and hence develop robust sensing algorithms. To partly 
counter these challenges, we develop a two-stage sub-space algorithm that detects multiple narrow-band 
analog frequency-modulated signals generated by WMs. The performance of the algorithm is verified by 
using experimentally captured low power WM signals with received power ranging from -100 to -105 
dBm. The problem of differentiating between the WM and other narrow-band signals is left as a future 
work. 
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I. Introduction 

Classical solution to the problem of spectrum allocation has been to divide it into non- 
overlapping frequency bands and assign each band to a wireless user/technology. While this 
avoids interference between users, it does not necessarily result in efficient spectrum utilization pi. 
This triggered a lot of interest in the wireless communications community to develop efficient 
methods that dynamically utilize the unused spectrum [|2|-||5j. In particular, the significant 
white space present in the TV bands (VHF/UHF) received a considerable attention and led to 
the conception of IEEE 802.22 wireless regional area networks (WRAN) ||6), [|7J. The basic 
underlying idea is to allow unlicensed users to use the licensed spectrum without causing harmful 
interference to the licensed users. One of the ways for the unlicensed users to achieve this is to 
identify the white spaces by spectrum sensing and then evacuate the spectrum when licensed user 
tries to access it. While spectrum sensing was a requirement for unlicensed TV Band Devices 
(TVBDs) in the FCC's Second Report and Order issued in 2008 [8], devices which access a TV 
Band Database are no longer required by law to perform spectrum sensing ||9). In stead, they are 
legally required to only perform geolocation in order to avoid causing harmful interference to 
licensed users of the spectrum. This is perhaps due to the limited performance of the sensing 
prototypes received by the FCC, and hence the need of efficient spectrum sensing still remains 
central to the success of decentralized DSA. 

The cognitive use of the TV Band requires protecting licensed ATSC, NTSC, and Broadcast 
Auxiliary Service (BAS) devices including Wireless Microphones (WMs). While there has been 
significant progress in developing detection methods for TV signals due to their defined structure, 



it has not been the case for WMs [10]-[12|. The main challenges has been the low transmit 
power of the WMs, lack of technical standard among various manufacturers and the resemblance 
of certain WM signals to narrow-band interference signals, such as spurious emissions fl"3] |. 
WMs typically operate with a transmit power ranging from 10 mW to 250 mW, though typical 
professional WMs operate with a transmit power of around 50 mW. Since WMs are used by 
live performers, their location is frequently in motion, and hence the field strength at any given 
point in space varies rapidly as a function of time due to fading and multipath interference. 
Almost all WMs use analog Frequency Modulation (FM), with a channel bandwidth of no greater 



than 200 kHz [14|. WMs may also operate on any carrier frequency that is a multiple of 25 
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kHz, and 25 kHz from the edge of the TV bands. Along with these uncertainties, WM signal 
detection is further complicated by real-world transmitter designs. Because of their relatively 
narrowband nature, nonlinear active devices such as amplifiers can generate Intermodulation 
Distortion (IMD). Such signals may be generated in WM transmitters (in which case, the IMD 
products are transmitted over the air), or in TVBD receivers (in which case the IMD products 
are only seen by the TVBD). In either case, the IMD products are nearly indistinguishable from 
the actual signals. The only way to mitigate this problem is to design receiver systems with a 
wide dynamic range to avoid operating in a nonlinear region. 

Under these uncertainties, it is extremely difficult to abstract features of WM signals and hence 



develop robust sensing algorithms that meet the suggested guidelines [15|, [16|. Matched filtering 
is optimal detection strategy if the TVBD receiver has a priori knowledge of the WM signal and 



the effective channel, which is not practical in general [17|. Energy detection is the simplest 
sensing algorithm to implement but it is not well suited for WM signals due to two main reasons: 
1) it employs non-coherent detection and hence requires longer sensing time to achieve the 
probability of detection constraints, and 2) the detection thresholds are highly susceptible to the 
changing noise levels. In fact, it is not possible to meet sensing requirements with a single energy 



detector under noise uncertainty, even if the sensing is done for a long time [ 18]— [20]. Another 
popular method is to estimate the spectral density of the received signal using maximum entropy 
principle [|2T|-[[23|. While this principle is analytically appealing, it has not been applied yet to 
WM sensing, perhaps due to the challenges in modeling the practical WM signals accurately. 
This is an interesting open problem but since modeling is not the main focus of this paper, we 
leave it as a possible future work. Recently, some interesting approaches for WM signal detection 



have been proposed in the literature, e.g., [26| proposes a cyclostationarity-based detector, (27 



proposes a Singular Value Decomposition (SVD) based detector, [28 1 proposes an Eigen value 
bases detector, and J29| proposes a covariance based detector. While all these methods perform 
well under their respective simulation models, their performance in real-world WM sensing is 
not clear. Moreover most of these methods in their current form are capable of detecting only a 
single WM signal and the extension of the analysis to multiple WM case does not seem to be 
trivial. 

Motivated by these recent developments, we propose a two-phase SVD based algorithm to detect 
multiple WM signals. In the training phase, we first compute the average noise autocorrelation 
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matrix from the known sets of noise data, which is then used to whiten the autocorrelation 
matrices of all the known sets of noise, unmodulated WM signals and modulated signals. These 
whitened autocorrelation matrices are then used to determine the empirical distributions of the 
test statistics, which are in turn used to set thresholds for signal detection. In the detection phase, 
we compute the autocorrelation matrix of the received signal and whiten it by subtracting the 
average noise autocorrelation matrix computed in the previous phase. The SVD of the resulting 
autocorrelation matrix is then computed to find the number of WM signals present. Based on 
the number of signals, the space is decomposed into signal subspace and noise subspace. The 
autocorrelation of the signals is recovered from the signal subspace, which is then used to find 
the center frequencies of the WM signals. We distinguish our work from other recent works, 



especially [27 1, [28 1, in two major ways. First, we extend previously proposed sub-space methods 
to the multiple WM sensing case by employing multiple test statistics. Second, we use real signals 
captured over-the-air to verify the detection performance of the algorithm. We use pass-band 
signals with 6 MHz bandwidth to closely emulate the television channel. The experimental data 
used for verification is representative of harsh channel conditions with received power ranging 
between -100 to -105 dBm. It is important to note that the problem of differentiating WM signals 
from other narrow-band signals is not considered in this work and may require techniques from 



patter recognition or robust statistics for an efficient solution [30|, [31 1. 

The paper is organized as follows. The noise characteristics are briefly studied by using 
over-the-air signals in Section II. As expected, the noise comes out to be stationary with a near- 
Gaussian distribution, but is not white. The proposed SVD based two-stage sensing algorithm is 
discussed in Section III. The detection performance of the proposed algorithm is numerically 
studied in Section IV. We first study the detection and false alarm probabilities using computer 
simulation and then verify the working of the algorithm on real over-the-air WM signal data. 
The paper is concluded in Section V. 

II. Measurement Campaign and Noise Characterization 

For the purpose of noise characterization and algorithm performance evaluation, different 
data sets were generated containing unmodulated WM signals, modulated WM signals, and 
background noise respectively. The modulated and unmodulated signal sets had a WM carrier at 
8 MHz, and were sampled at 33 MHz. The signal was bandpass-filtered to simulate a 6 MHz 
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television channel. The received power level in these data sets was specified as being between 
—100 dBm and —105 dBm. The measurements were captured inside a building and contain a 
WM transmitter located at varying distances from the receiver. Not only does this provide a 
very wide range of received signal strength signals, but it also accurately replicates real-world 
conditions such as multipath and external interference. 

We use three noise data sets, referred henceforth as NoiseData{l:3}, to examine various 
properties of the background noise. The histogram of the noise samples, shown in Fig. [TJ 
appears to be a truncated Gaussian, although we did not formally test it for Gaussianity. With 
the assumption that the noise is Gaussian, we can simplify the complexity of the algorithm 
significantly. The variance of the noise for NoiseData{l:3} is shown in Fig. [2j The plot indicates 
that the noise variance does not vary significantly as a function of time. A similar observation 
was made with the noise mean. The observations are consistent across the noise data sets. 

Figs. [3] and [4] show the autocorrelation and power spectral density (PSD) of the noise samples. 
The PSD clearly shows that the noise was band-pass filtered prior to sampling. The autocorrelation 
plot reflects this as well, showing correlation up to approximately 15 samples. In addition, the 
noise correlation properties do not vary significantly across the noise data sets. 

Figs. [5] and [6] show the autocorrelation and PSD of the data samples containing a modulated 
WM signal. As expected, the autocorrelation differs significantly from Fig. [3] in that the correlation 
with signal present extends out far past 15 samples. It is this property that allows us to use a 
correlation-based technique to determine whether a WM signal is present in a given data set or 
not. The autocorrelation of the data samples containing a silent WM signal resembles that in 
Figure [5J with the exception that the correlation is lesser. 

These observations are in line with the intuition and we conclude that the noise is stationary 
with a near-Gaussian distribution, but is not white. Rather, the noise is band-limited (colored) 
due to filtering and sampling in the receiver. 

III. SVD Based WM Sensing Algorithm 

The working of our proposed algorithm can be classified into two main phases, viz., training 
phase and detection phase. Main steps involved in these two phases are as follows: 
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Training Phase: 

1) Characterize the noise by using the training data sets NoiseData{l:3} containing no WM 
signal. 

2) Compute the average noise autocorrelation matrix from all the training data sets. 

3) Subtract this average noise autocorrelation matrix from the autocorrelation matrix of each 
NoiseData to get the whitened noise correlation matrices. 

4) Determine the thresholds of the test statistics by using training sets of noise and WM 
signal. 

Detection Phase: 

1) Compute the autocorrelation matrix of the received signal. 

2) Whiten the noise by subtracting the average noise autocorrelation matrix from the received 
signal correlation matrix. 

3) Find the SVD of the resulting correlation matrix and determine the number of signals 
present by comparing the test statistics with the thresholds. 

4) Decompose the space into signal sub-space and noise sub-space according to the number 
of WM signal present. 

5) Recover the autocorrelation of the signal from the signal autocorrelation matrix (signal 
sub-space). 

6) Find the center frequencies by the Fourier analysis of the signal autocorrelation. 

With this brief introduction to the proposed algorithm, we now discuss it in detail. The received 
signal samples are modeled as: 

x(n) = s(n) + r](n), (1) 

where s(n) represents the transmit signal and rj(n) models the cumulative effect of thermal noise 
and interference, referred henceforth as only noise. For a received signal sequence of length L, the 
Lx L autocorrelation matrix can be represented as R^. = R s + R T; , where R^, = E {x(ri) x H (n)} 
is the autocorrelation matrix of the received samples, R s = E {s(n) s H (n)} is the autocorrelation 
matrix of the transmitted signal samples, and R^ = E {r)(n) r] H {n)} is the autocorrelation matrix 
of the noise samples. Note that s H denotes the conjugate transpose of vector s. In the case when 
the noise is white, R,, = II, where II is an identity matrix of order L. The value of L is 
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chosen as a trade-off between complexity and performance. We will comment more on the choice 
of L in the numerical results section. 



A. Noise Whitening 

To perform noise whitening, we first compute an estimate of the average noise autocorrelation 
matrix R T; by taking average of the autocorrelation matrices of the three noise sets NoiseData{l:3}. 
This average autocorrelation matrix is then used to cancel off the noise correlation from the 
received signal autocorrelation matrix as follows: 

R s = Rr — R»; = Rs + R^t, (2) 

where R^ is the correlation matrix of the residual noise. R s is the estimate of the correlation 
matrix of the WM signal assuming that the correlation matrix R M of the residual noise is negligible 
compared to R s . This basic method should not be confused with the popular noise whitening 
methods where the correlated noise is passed through a filter to achieve zero mean and identity 
covariance matrix. There are two main reasons for choosing this basic method, despite it being 
not as accurate as the standard whitening method: 1) it is difficult to design the whitening matrix 
with the limited noise data sets, and 2) more importantly, it serves the purpose of removing the 
prominent interference signals that might be present in the ambience. 

B. SVD and Test Statistic Threshold 

The SVD method can be applied to the cleaned autocorrelation matrix of received signal R s 
to obtain 



R S = USY H = [U S U M ] 



s s 



o s M 



[V S V»} H , (3) 



where S s and S M are diagonal matrices whose values correspond to the singular values in the 
signal subspace and noise subspace, respectively. It is straightforward to show that each WM 
signal produces two non-zero singular values. Therefore, for N s WM signals, diag(S s ) will 
be [Ai, A 2 , . . . A 2 /v s ] and diag(S^) will be [A 2 /v s +i, A 2 at 3+2 , • • • Al], where the singular values are 
arranged in the decreasing order. E.g., if there is only one WM signal present, diag(S s ) = [A 1; A 2 ] 
and diag(S u ) = [A 3 , A 4 . . . X L ]. It should be noted that the singular values corresponding to 
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noise sub- space are all of same order and much smaller than those of the signal sub- space. The 
natural test statistic, therefore, is to take the ratio of alternate singular values, i.e., A2X-1/A2X+1, 
1 < x < (L — 2)/2. In case only white noise is present, we get A2K-1/A2X+1 — h V x. However, if 
N s number of signals are present, we get X^x-i/^x+i >> 1 for 1 < x < iV,, and A2X-1/A22+1 = 1 
for N s + 1 < x < (L - 2)/2. 

The threshold values, \ n for 1 < i < (L — 2)/2, are determined by applying SVD to the 
whitened autocorrelation matrices of known sets of noise, unmodulated WM signals and modulated 
WM signals. The empirical distributions of test statistics are determined in each of these cases 
and thresholds are set to achieve predetermined probability of false alarm (P/ a ) or probability of 
detection (P d ). Those familiar with the problem of determining empirical distributions are likely 
to immediately recognize that it is very hard to get a clean estimate from a very limited set of 
observations. Therefore, we identify the duration and frequency of training as two main practical 
parameters that affect the performance of the proposed algorithm. 

C. Center Frequencies of WM Signals 

After identifying the number of WM signals present, the space is decomposed into noise and 
signal subspace as shown in equation ([3]). The autocorrelation of the cleaned signal is recovered, 
which is nothing but the first row of the autocorrelation matrix R s . PSD is then determined by 
taking the Fourier transform of the signal autocorrelation. The locations of the WM signals can 
be estimated by the locations of N s highest peaks in the PSD. Since the autocorrelation matrix 
R s is reconstructed only from the signal subspace, the resulting PSD is much cleaner than that 
obtained from correlation matrix R x . or even R s . This makes determining the location of the 
WM signals much easier and more robust. 

IV. Numerical Results 

In this section, we determine the detection performance first by using computer simulation 
and then by applying the algorithm on real over-the-air WM signals. The computer simulation is 
used to study finer performance details, which is not always possible to do in case of a limited 
set of real-world data. 
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A. Computer Simulation 



We generate WM signals based on the procedure explained in [32]. Two different data sets 
are generated: the training set and the test set. The two data sets essentially consist of same type 
of data. The training set is used to obtain the thresholds for the test statistics, while the test set 
is used to measure the performance of the algorithm. The sampling rate used for the simulation 
is F s = 33.33 MHz and each test set is 20000 samples long. 

1 ) Generation of Training and Test Set: We generate both the noise-only data and the WM 
signals embedded in noise for both the training and the test sets. The noise signal is generated 
by filtering additive white gaussian noise (AWGN) with a bandpass filter of 6 MHz bandwidth 
centered at the IF of 8 MHz. Fig. [7] presents the PSD of noise-only signal. For data with WM 
signals, the number of signals is randomly varied from one to five. The locations of the signals 
were fixed to pre-specified frequencies of 6,7,8,9 and 10 MHz for ease of verification. For 
example, if the number of signals was selected to be two, the signals would be fixed at the 
frequencies of 6 and 7 MHz. If five were present, they would be centered at 6, 7, 8, 9 and 10 
MHz. The SNR of the WM signal is varied from -30 dB to -15 dB. For this paper, only the WM 
signal with loud speaker is simulated. The silence and soft speaker modes are left as a future 
work. Fig. |] shows the PSD of WM signals at 6 and 7 MHz at SNR of -20dB. Fig. [9] shows the 
PSD of wireless microphone signal with loud speaker that was used for the simulation. 

2) Test Statistic Threshold Determination: In order to derive the threshold for the test statistics, 
the distribution of the test statistics was first obtained. As shown in Fig. [TO], at relatively high 



SNR of -15dB, the distributions of test statistics for noise-only and signal plus noise cases are 



well separated. However, Fig. 11 shows that at relatively low SNR of -25 dB, the distribution 
of test statistics of the two cases are not well separated, which is inline with the intuition. The 
thresholds are determined using simple hypothesis testing. 

H , No wireless microphone present 
Hi, Wireless microphones present 

Given the distributions of the test statistics, we find A r such that Pf a = P[X > X T \H ] < 0.1, A is 
the test statistic, and A T is the threshold. However, since the distributions do not perfectly match 
the true probability distribution of the test statistics, the threshold was later adjusted iteratively 
by finding the actual Pj a . Five thresholds were determined for each SNR value ranging from -30 
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dB to -15 dB. Table U shows the thresholds obtained for each SNR. 



3) Detection Results: Fig. 12 shows the probability of detection and probability of false alarm 
based on the derived thresholds for AWGN case. The detection results for the Rayleigh fading 
case will be included in the future version of this paper. The results show that the algorithm 
successfully detects multiple wireless microphone signals with high probability at SNR as low as 
-23 dB if a sufficiently long signal sequence is taken. Although not perfectly, the probability of 



false alarm is also approximately bounded by 0.1. Fig. 13 shows the PSD of the original signal 
at SNR of -25dB and the PSD obtained after processing the original signal. It shows that five 
signals at 6, 7, 8, 9 and 10 MHz are correctly detected. 

4) Effect of Varying L: Calculating the autocorrelation matrix, which is one of the steps 
required in the algorithm, is computationally expensive for longer sequences (higher values of 
L). Furthermore, working with large matrices is typically not desirable for real-time sensing 
algorithms. Throughout the paper, we have considered the value of L to be 500, which was 
chosen to make sure that the resulting signal sequence was long enough to achieve the best 
possible detection performance. We now vary the value of L and study the trade-off between 



performance and complexity in Fig. 12 As expected, the decreasing value of L degrades the 
detection performance of the algorithm. We note that L ps 200 is a good tradeoff between 
performance and complexity. 

B. Using Real Data 

The proposed algorithm has been verified and tested on real over-the-air data. The empirical 
thresholds that were used for detection are as follows: X n = 1.8741, A T2 = 1.4505, A T3 = 



1.5743, A r4 = 1.4806 and A T5 = 1.3152. Fig. 14 shows the PSD obtained from the cleaned-up 



autocorrelation function R s for a data set containing WM signal. As mentioned earlier, it is not 
practical in this case to estimate detection and false alarm probabilities due to the limited data 
set. Therefore, we analyze the data sets by visual inspection as an extra step independent of the 
proposed algorithm. Interestingly, the visual inspection has revealed that a few data sets contain 
more number of signals than are detected by the proposed algorithm. Technically this means that 
some of these signals are not leading to independent columns of the autocorrelation matrix, and 
hence the number of non-zero singular values is lower than expected. One possibility is that there 
could be intermodulation products present in the received signal that share the same subspace as 
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the WM signals itself, in which case the SVD algorithm may not observe any singular values 
corresponding to the intermodulation products. 

V. Conclusion and Future Direction 

Detecting WM signals in poor SNR conditions using test equipment is a difficult task. Doing 
the same task with a consumer-grade receiver over a wide band will prove to be even more 
challenging. This paper has shown how a simple sub-space method can detect the presence of 
extremely weak WM signals and determine what frequencies they operate on. 

The current version of the proposed algorithm faces one drawback. It is not capable of 
examining each signal individually to determine whether it exhibits the characteristics of a 
frequency-modulated analog audio signal. The next step in this process is to incorporate this 
information into the algorithm, as well as the capability to calculate potential IM3 products and 
compare the results to the detected signals to reduce the possibility of false positives (especially 
over multiple TV channels). More measurement campaigns have to be conducted under different 
wireless environments in order to obtain more statistically accurate calculations of the test statistic 
thresholds. One other important research direction would be to evaluate the complexity and 
performance of the proposed algorithm on software radio platforms including small form factor 
handsets. 
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TABLE I 
Test Statistic Thresholds 



SNR (dB) 




A T 2 


A T 3 


A T 4 


A T 5 


-30 


1.568 


1.252 


1.251 


1.196 


1.143 


-29 


1.655 


1.294 


1.221 


1.194 


1.174 


-28 


1.505 


1.357 


1.205 


1.202 


1.168 


-27 


1.657 


1.290 


1.236 


1.198 


1.205 


-26 


1.459 


1.304 


1.260 


1.166 


1.217 


-25 


1.569 


1.288 


1.242 


1.228 


1.157 


-24 


1.463 


1.387 


1.229 


1.231 


1.227 


-23 


1.678 


1.319 


1.220 


1.205 


1.176 


-22 


1.517 
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Fig. 1. Histogram of background noise computed from measurement data. 




Fig. 2. Variance of background noise at different time instances. 
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Fig. 3. Autocorrelation of band-limited background noise samples. 




Frequency (MHz) 



Fig. 4. PSD of band-limited background noise. 
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Fig. 5. Autocorrelation of loud (modulated) WM signal with colored noise. 
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Fig. 6. PSD of modulated WM signal with colored noise. 
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Fig. 9. PSD of simulated wireless microphone signal: loud speaker. 
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Fig. 10. Histogram of test statistic at SNR = 



-15 dB. 
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Fig. 11. Histogram of test statistic at SNR = -25 dB. 
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Fig. 12. Detection result for simulated data in AWGN. 
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Fig. 13. PSDs of the original data and the cleaned-up autocorrelation: five signals present. 




-180 - 

-200 I '< '< '< ' ' ' 1 

2 4 6 8 10 12 14 16 

Frequency (MHz) 



Fig. 14. PSDs obtained from the cleaned-up autocorrelation function R s . 



